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Abstract 

An algorithm to compute a good basis of the Brieskorn lattice of a cohomologically 
tame polynomial is described. This algorithm is based on the results of C. Sabbah 
and generalizes the algorithm by A. Douai for convenient Newton non-degenerate 
polynomials. 
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Introduction 



Let / : C""*"^ > C with n > 1 be a coliomologically tame polynomial func- 
tion [1]. This means that no modification of the topology of the fibres of / 
comes from infinity. In particular, the set of critical points C(/) of / is finite. 
Then the reduced cohomology of the fibre f~^{t) for t ^ C(/) is concentrated 
in dimension n and equals where fi is the Milnor number of /. Moreover, 
the n-th cohomology of the fibres of / forms a local system if" on C\D(/) 
where D(/) = /(C(/)) is the discriminant of /. Hence, there is a monodromy 
action of the fundamental group ni(C\D(/),t) on if". 

The Gauss-Manin system M of / is a regular holonomic module over the Weyl 
algebra <C[t]{dt) with associated local system ii" on C\D(/). The Fourier 
transform G := M of M is the C[r](9T-)-module defined by r := dt and 
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dr — —t. The monodromy Too of M around D(/) can be identified with the 
inverse of the monodromy Tq of G at 0. It turns out that dt is invertible on 
M and hence G is a C[t, 6']-module where 6 := A finite C[t]- resp. C[^]- 
submodule L C G such that L[9] — G resp. L[r] = G is called a C[r]- resp. 
C[^]-lattice. The regularity of M at oo imphes that G is singular at most in 
{0, oo} and where := {r = 0} is regular and oo := = 0} of type 1. In 
particular, the V-filtration V, on G at consists of C[t] -lattices. 

The Brieskorn lattice Gq C G is a ^-invariant C[^]-submodule of G such that 

G = Go[r]. C. Sabbah [1] proved that Gq is a free C[t]- and C[^]-module of 
rank fi. In particular, G is a free C[r, ^^]-module of rank /v,. By definition, the 
spectrum of a C[^^]-lattice L C G is the spectrum of the induced V-filtration 
V,{L/9L) and the spectrum of / is the spectrum of Gq. 

C. Sabbah [1] showed that there is a natural mixed Hodge structure on the 
moderate nearby cycles of G with Hodge filtration induced by Go- This leads to 
the existence of good bases of the Brieskorn lattice. For a basis = 0i, . . . , 0^ 
of a ^-invariant Cf6'l-lattice, 



where A± e €[e]^''"'. A C[^]-basis of Go is called good if A± ^ + OAj 
where A^AjeC'"'''', 



and 0j e V^jGo for all i E One can read off the monodromy = Tq^ 

from A!t immediately. The diagonal a = ai, . . . , is the spectrum of / and 
determines with gi^ Aq the spectral pairs of /. The latter correspond to the 
Hodge numbers of the above mixed Hodge structure. 

Analogous results to those above were first obtained in a local situation where 

/ • (C", 0) > (C, 0) is a holomorphic function germ with an isolated critical 

point [2,3,4,5,6,7]. In this situation, the role of the Fourier transform is played 
by microlocalization and the algorithms in [8,9] compute Aq and Ai for a 
good C{{6'}}-basis of the (local) Brieskorn lattice. But [9] and [8, 7.4-5] do 
not apply to the global situation. 

A. Douai [10] explained how to compute a good basis of Gq if / is convenient 
and Newton non-degenerate using the equality of the V- and Newton filtration 
[11,1] and a division algorithm with respect to the Newton filtration [12,13]. 

The intention of this article is to describe an explicit algorithm to compute 
a good basis of Go for an arbitrary cohomologically tame polynomial /. This 
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algorithm is based on the following idea: 

Let coordinate system on C""*"^. Then the Brieskorn lattice 

Go can be identified with the quotient 

/ " 

€[x,0]/Y.{d,,{f)-0d,,)i€[x,e]) 

' 1=0 

of non-finite C[^]-modules. The degree with respect to x defines an increasing 
filtration €[x, 9], by finite C[^]-modules on €[x, 9] and hence 

G'o'' C[x, 9], / (C[x, 9]^ n - Od^Mx, 9]i)) 

are finite C[^]-modules. For A; » and / ^ 0, Gg ' — by the finiteness of 
Go- But, a priori, there is no bound for these indices. 

By Grobner basis methods, one can compute cyclic generators of a 
invariant C[0]-sublattice Gq'' C Gq. By an argument of A. Khovanskii and 
A. Varchenko [11], Gq'' = Go if and only if the mean values of the spectra 
coincide. By the t-invariance of GS'' , one can compute the spectrum of Gg' 
like that of Gq below. The mean value of the spectrum of Gq is known to be 
So if the mean value of the spectrum of Gg'' is not then one has to 
increase k. This process terminates with Gg'' = Gq. 

Then one can compute for the C[^^]-basis of Gg. By a saturation process, 
one can compute the V-filtration and, by a Grobner basis computation, the 
spectrum of Gq and the Hodge filtration. Then one can compute a C[t, ^]- 
basis of G which is compatible with the V-filtration refined by an opposite 
Hodge filtration. In terms of this basis, one can compute a good basis of Gg 
by a simultaneous normal form computation and basis transformation. 

We denote rows vectors v by a lower bar and column vectors v by an upper 
bar. In general, lower indices are column indices and upper indices are row 
indices. We denote by {M} the set and by {M)R the i?-linear span of the of 
columns of a matrix M . We denote by lead the leading term and by lexp the 
leading exponent with respect to a monomial ordering. We denote by E the 
unit matrix and by Cj the iih. unit vector. 



1 Gauss-Manin system 

Let / : C**"*"^ > C with n > 1 be a polynomial function. Let O be the sheaf 

of regular functions and (fi', d) the complex of polynomial differential forms on 
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(|-in+i 'jj^gj^ ii^Q Gauss-Manin System /+0 of / is represented by the complex 
of left C[i](at)-modules 

(Q'+"+i[5,],d-5id/) 

[5, 15] and has regular holonomic cohomology [14, VII. 12. 2]. The coefficients 
of the differentials are the differentials of the complexes (i^*,d) and {fl*,df). 

Lemma 1 (Poinceire Lemma) The complex of C-vector spaces 

>C >(n%d) ^0 

is exact [15, Ex. 16.15]. 

From now on, we assume that set of critical points C(/) of / is finite. Then 
the following lemma holds. 

Lemma 2 (De Rham Lemma) 

(1) H*^(f]*,d/) =0 for k + 

(2) dimcH"+^(0',d/) < oo. 

Proof. If C(/) is finite then 

€[x]/{d{f)) ^ Q"+7d/ A Q"-^ = H"+^(Q',d/) 

is a finite C-vector space and hence d{f) is a regular sequence in <C[x\. Then 
the cohomology of the Koszul complex {Q*,df) is concentrated in dimension 
n + 1 [15, Cor. 17.5]. 



The image Mq of Q""*"^ in M :— H°(/+0) is the key for an algorithmic approach 
to the Gauss-Manin system. It determines the differential structure of M and 
it can be identified with a quotient of C[x]. 

Proposition 3 

(1) H'=(/+O) = 0/orA;^{-n,0}. 

(2) dt is invertible on M. 

(3) Mo = 1]"+Vd/ Adl]"-i. 



Proof. This follows from Lemma 1 and 2 and [5, 15.2.2]. 
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The Fourier transform M of M is the left C[T](9T-)-module defined by the 
isomorphism 

r := dt, dr := -t 

of €[T]{dr) and €[t\{dt) [16, 2.1]. By Proposition 3, M is a C[e](ae)-module 
with 

e ■= T-\ de := -T^dr 
and Mq is a C[6']-submodule. Note that t = 

Definition 4 Let G be the €[9]{d0) -module M. Then the Brieskorn lattice Go 
of f is the C[9]{t)-submodule Mq ofG. 

Since M is regular at oo, G is singular at most in {0, oo} where := {r = 0} 
is regular and oo := = 0} of type 1 [17, V.2.a]. 

From now on, we assume that / is cohomologically tame. By definition [1, 8], 
this means that there is a compactification 




where C"^+^ is quasi-projective and / is proper such that, for all t G C, the 
support of the vanishing cycle complex (/>j_jRj*Q is a finite subset of C"'"'""^. 
In particular, C(/) is finite and hence, by Lemma 2, the Milnor number 

II := dime H"+^(n', d/) = dimc(l^"+Vd/ A Q") 

of / is finite. Then the following theorem holds. 

Theorem 5 (C. Sabbah [1, 10.1-3]) Go is a free €[t]- and C[e]-module 
of rank /i. 

In particular, G is a free C[t, ^]-module of rank /j,. 



2 Brieskorn lattice 

A finite C[t]- resp. C[6']-submodule L cG such that L[9] = G resp. L[r] = G 
is called a C[r]- resp. C[^]-lattice. By Theorem 5, a lattice is free of rank fi. In 
terms of a C[^]-basis of Gq, the C[^^](9e)-module structure of G is determined 
by the basis representation of t on Gq. The following lemma shows that the 
latter is determined by a matrix with coefficients in C[^]. 
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Definition 6 Let (p be a basis of a t-invariant (C[6]-sublattice L <Z G. Then 
the matrix e C[^]''^'* oft with respect to is defined by 

^At := t^. 

Lemma 7 Let (f) be a basis of a t-invariant €>[9]-sublattice L <Z G. Then 
Proof. Since t = O'^de, 

toJj2p,e')=tj:iP,e' 

^ k ^ k 

^J2t{^p,)e' + ^p,e'doe' 

k 

^lo{A± + 9%)(^^p,9') 
and hence i o = o {A^ + O'^de) . 

The following lemma gives a presentation of the C[^](i)-module Gq. This 
presentation shall be used to compute t on Gq. 

Lemma 8 There is an isomorphism of C[9]{t) -modules 

Go^n^^'[e]/{df-ed){n''[e]). 

Proof. By Theorem 5, 

G = q''+^[t, e]/{df - ed){Q''[T, e]). 

By definition, Gq is the image of in G and hence 

Go = Q''+^[e]/({df - M)(0"[r, e]) n n^'+^ie]). 

By Lemma 2, dker(d/) C d/ A dQ"""^ C ker(d/) and hence 

(d/ - ed){Q''[T, e]) n ^""+^[6] = (d/ - ed){Q''[e]). 

Let be coordinates on C"'*"-'^ with corresponding partial deriva- 

tives 9 = da;o, ■ ■ ■ ,dx„. Let 

t:^f + e%e€[x,e]{de). 
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Then, by Lemma 8, we can identify 

G = c[x, T, e]/{d{f) - ed){c[x, T, er+') 

as C[t, ^](i)-modules and 

Go^(C[x,e]/{d{f)-ed){€[x,er+') 

as C[^^](t)-modules. These modules are quotients of non-finite C[^]-modules. 
On the numerator and denominator, the degree with respect to x defines an 
increasing filtration by finite C[^]-modules. The following algorithm computes 
t on Go by an approximation process with respect to these filtrations. 

Definition 9 The degree deg^ with respect to x defines an increasing filtration 
€>[x,9], on€>[x,9] by finite C[9] -modules 

C[x,e]k :^{peC[x,9]\deg^{p)<k} 

such thatt€[x,9], C €<[x, 9],^deg{f) ■ define the finite C[9]-raodules 

G'o := C[x, 9]k/{{d{f) - 9^)(C[x, 9r+') n C[x, 9]^), 
G'o' := €[x, 9]k/l{d{f) - mi^k, ^ir') n C[^, 9]^). 

Algorithm 1 

Input: (a) A cohomologically tame polynomial / G C[x]. 

(h) An integer k >0. 
Output: (a) A vector G C[x, such that [(/)] is a basis of a t-invariant G[9]- 
lattice Lk C Go o,nd Lk = Go for k ^ 0. 
(b) The matrix A = A^^ e G[9f''^'. 

(1) Setl:^k. 

(2) Set I -.= 1 + 1. 

(3) Compute a reduced Grobner basis 

g GB({d{f) - d9){jf-^ei) \ e x [0,n], |a| < /) 

of (d(f) — d9) (Clx, 9]^'^^^ with respect to a monomial ordering > on 
{x^' I (a, i) e ]N"+^ X IN} such that 

\a\ > 1^1 =^ a > ^, 
{a,i) > <^ a> PV {a = PAi > j) 

for all {a,i), {f3,j) e lN"+i x IN. 

(4) Find the minimal ko with 

ko<\a\<k^x^e (lead(£))C[e]. 
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(5) Compute (f) e such that [0] are cyclic generators of 

Co'"'' = €[x,e]J{g, I deg,(^,) < ^o)C[^] = Go''' 

and p := rk(G'o"'') using [18, 2.6.3]. 

(6) Ifp>iior^>p — II then go to (2). 

(7) If p < 11 then set k := k + 1 and go to (2). 

(8) Ifko + dcg(/) > A; then set k := k + 1 and go to (2). 

(9) If [0] is not a C-basis of C[a;]/(9(/))C[a;] then set k := k + 1 and go to 
(2)7 

(10) Compute a normal form NF{t(j),g) of tcf) with respect to g. 

(11) Compute the basis representation A e C[^]'^^'* o/ [NF(i0, gi)] with respect 
to the C[9]-basis [0] of Lk := Gq°'^. 

(12) Return 4> and A. 

Lemma 10 Algorithm 1 terminates and is correct. 

Proof. Since d{f) — d9 is C[6']-linear, 

{{d{f) - d9){x^ei) I {a,i) e 1N"+^ x [0,n], \a\ < l)C[9]. 
By definition of the monomial ordering, 

im - mCCk, 9]^') n C[x, 9], = {g, I deg,(^,) < k)€[9] 

and hence, by definition of kg, 

G'^''^€[x,9y{g^\deg^{g,)<k)€[9] 

= <C[x,9U/{g, I deg,(^,) < ko)C[9] = G^'^ 

Because of step (2), I is strictly increasing for fixed k. There are C[^]-linear 
maps 

Gl^'^G^^^Go 

where Lk is an isomorphism for /c 3> and vr^ ; is an isomorphism for fixed 
k and / ^ 0. By Theorem 5, Gq is a free C[6']-module of rank at most p. 
Hence, if condition (6) holds then tt*^'' is not an isomorphism and if condition 
(7) holds then i*^ is not an isomorphism. 

By Theorem 5, there is a '0 G C[a;, ^^]^ such that [ip] is a C[^]-basis of Gq. In 
particular, is an isomorphism for k > degj.{ilj) and hence 
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is an isomorphism and conditions (6) and (7) do not hold for k > deg^ (■?/') and 
/ > 0. By Lemma 8, for each a e IN"'^^ there is a matrix M- e C[6'p^'' such 
that 

- e id{f ) - d9){C[x, 9]''+^). 
If |a| > deg^(^) then x^-i^M^ e id{f) - de){C[x, 9]"^) D €[x, and hence 

e lead((9(/) - d9){<Ck: ^1^')) = {lee.d{g))C[9] 

for / ^ 0. Hence, by definition of ko, ko < deg3.(V^) for k > deg^{ip) and / 3> 
and, in particular, condition (8) does not hold for k > deg^{ip) + deg(/) and 
I » 0. Since [ijj] is a C[^]-basis of Go, [V^] is a C-basis of 

Go/9Go = C[x]/{d{f))€[x] 

and hence condition (9) does not hold for k > deg^(^) and I » 0. This proves 
that the algorithm terminates. 

Since [0] is a C-basis of €[x}/ {d{f))€[x\ = Go/9Go, o tt'"''' is injective and 
[0] is a basis of the C[^]-lattice Lk = Gq"'' C Gq. Since e C[ai,^]^^ and 
ko + deg(/) < k, t(j)_e C[x, 6*]^ and hence, by definition of ko, NF{t^,g) e 
C [x, 9] 1^ . By Lemma 8, 

m = M = [NF(t0,5)] = [M] = [l]A 

and hence Lk is ^-invariant and A — A^^. This proves that the algorithm is 
correct. 



A priori, we do not know a ko such that = Go for all k > ko. We shall solve 
this problem by a criterion on the spectrum with respect to the V-filtration. 



3 V-filtration 

Definition 11 The V-filtration V, on €[T]{dr) is the increasing filtration by 
VoC[t] {Ot-) -modules 

Vk+lC[r]{dr) := VkC[r]{dr) + drVk€[r]{dr) 

for all k>0. 

Proposition 12 There is a unique V,<C[T]{dT-) -good filtration V, on G by 
€>[T]-lattices such that rdr + a is nilpotent on gr^ G for all a. 



9 



Proof. Since G is regular at 0, this follows from [19, 2.3.2, 4.1, 5.1.5]. 

Definition 13 V,G is called the V-filtration on G. 

The following criterion shall be used to compute the V-filtration on G. 

Lemma 14 Let L <Z G be a rdr -invariant €[T]-lattice with 

spec{-Tdr e End(L/TL)) (z[a,a - 1) 
for some a. Then L — V^G. 

Proof. Let spec (—tc^t- e End(L/rL)) = {a} with 

Let 4> : L/tL ¥ L be a C[r]-basis of L and 

:= (/)(ker((T9, + a^f e End(L/TL))) 

for all i e [1,1^]. Let 

u 
i=j 

for all i e [1,1^] and p e Z. Then [/, is an increasing filtration on G by tOt-- 
invariant C[r]-latticcs. By construction, rdr + ttj — p is nilpotent on gr^^^G 
and Uai-p — T^Ua^ for alH e [1, i/] and p e Z. Since 

= T^'-\Tdr+p- 1) 0Ca, © T^'lra, +p)L 

C T^'-l (ra, + P - 1) © T^irdr + p) + 
i=j i=l 

j-1 

Uaj-p+1 — ^ Cq,. © Cq. © T^"*""^!/ 

i=j i=l 
V j-1 

© ''"^ + Uaj-p, 
i=j i=l 

-p+i for p > ttj + 1 and hence U, is V,€[r]{dr)-good. 
Then, by Proposition 12, U,G = V,G and hence L = V^G. 

The following algorithm computes the V-filtration using the criterion in Lemma 
14. For a given C[6']-lattice with C[6']-basis (p, L := (0)C[r] is a C[r]-lattice 
with C[r]-basis and -rdr^ = where B = tA^ e Cir,^]/^^^. By a sat- 
uration process of L with respect to rdr, L is replaced by a r^^ -invariant 
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C[T]-lattice and is modified such that B G C[r]''^^. Then a sequence of 
basis transformations modifies such that spec(i^o) C [a,a — 1) for some a. 

Algorithm 2 

Input: The matrix A — e C[^]'*^'' for a basis (p of a t-invariant €[9]-lattice 



Output: (a) A matrix U e C[6']^^'^ such that 4>U is a C[T]-basis of Va for some 



(b) A matrix B = Ei>oBiT' E C[r]'*^'' such that -rdr{^U) = ^UB and 
spec(-Bo) = {«} with a>ai>--->a^>a — 1. 

(1) (a) Set k:=0 and Uq := E E C^^^. 

(b) Until {{tA - Tdr)iUk)} C {Uk)€[T] do: 

(i) Set k:=k+l. 

(ii) Compute Uk E C[^]''^^ with deg{Uk) < A;(deg(A) - 1) such that 

{Uk+l)C[T] = {Uk)C[T] + {{tA - Tdr){Uk))€[T]. 

(c) Set U := Uk- 

(2) (a) Set B = E.>o B,r^ := U'\tA - rdr){U) E C[t]''X'^. 
(b) Compute {a} := spec(i?o) and j E such that 



LcG. 



a. 



ai > ■ ■ • > aj > ai — 1 > aj+i > • ■ ■ > a^,. 



(c) If j = u then return U and B. 

(d) Compute Uq E GL^(C) such that 




where 



spec(5o' ) 
spec(Bo'^) 



{«!, . . . ,aj}, 




0. 
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(f) SetU := (^Ui T-^U2^ and 



B 



(g) Set ai :^ ai + 1 for i ^ j + 1, . . . ,1/. 

(h) Reorder a and redefine j e such that 



CKi > • • • > CKj > tti — 1 > Ctj+l > • • • > ai,. 



(i) Go to (2c). 
Remark 15 

(1) IfA^Ao + Ml then Uk ^ Uq ^ E. 

(2) If 

B = 




with spec(So' ) = {ai, 
can choose 



, aj} and spec(So'^) = {oij+i, ■ ■ ■ , Oii,} then one 
Uo 



E 



Lemma 16 Algorithm 2 terminates and is correct. 



Proof. 

(1) By Lemma 7, 

(0t/fc+i)C[T] = {m)€[T] + (0O (rA - Tdr){Uu))€[T] 
= (0f/,.)C[r] + r(0 o (A + e^de){Uk))€[T] 

= {±Uk)€[T]+Tdr{±Uk)€[T] 

and hence {(0f/fc)C[r]}fc>o is an increasing sequence of finite C[r]-modules. 
Since (0f/o)C[r] = €[tY , one can choose Uk E C[^]''''^. The V-filtration 
on G consists of finite and hence Noetherian re?,- -invariant C[T]-modules. 
For some a, {0} C V^G and hence (0[/fc)C[r] C Va for all > 0. This im- 
phes that the sequence {(t/fc)C[r]}fc>o is stationary. Then {^U)€[t] C G 
is a rc^T-invariant C[r]-lattice. 

(2) By Lemma 7, —tOt ocf) — rtocj) — (j)o {jA — tOt) and hence 

-rdr o^U^lUo(B- Tdr). 
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The T9T--invariance of the C[t] -lattice (0C/)C[t] is preserved since 



E 



V 



( ^1,1 ^-1^1,2^ 











(e 


\ 


1 

























and B^''^T~^ e C[tY^'^^~^\ The index j is strictly increasing since 



spec(i?o) = spec 



By - E r-'B^^ 

2,2 

/ 



B, 



= {«!, . . . , aj, aj+i + 1, . . . , a,, + 1} 

and hence the algorithm terminates. Then L := (0f/)C[r] C G is a rdr- 
invariant C[T]-lattice with spec(— rc?^ G End(L/rL)) C [a,a — 1) for 
a :— ai. Hence, by Lemma 14, L — Va and (j)U is a C[r]-basis of VaG. 



4 Spectrum 

The spectrum with respect to the V-filtration shall be used to check equahty 
of C[^]-lattices. 

Definition 17 

(1) The spectrum spec(F,) : Q >-]N of an increasing filtration F, on a 

finite vector space V is defined by 

spec(F,)(Q;) := dim(gr^* V) 

for all a E Q,. The spectrum spec(F,) of a decreasing filtration F* on V 

is defined analogously. 

(2) The spectrum of a G[6]-lattice L G G is defined by 

spec(L) := spec{V,(L/eL)). 

(3) The spectrum of f is defined by 

spec(/) := spec(Go)- 

The following algorithm computes the spectrum of a t-invariant C[^]-lattice 
by computing a Grobner basis compatible with the V-filtration. 
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Algorithm 3 

Input: (a) A matrix B = Y.i>oBiT'- G C[r]''^^ such that —Tdr(p = 4>B for a 
G[T]-hasis (j) of Va and spec(i?o) = {«} with a > ai > • • • > ai, > 
a — 1. 

(b) A matrix M G €[t,9]^^^ such that 0M is a basis of a t-invariant 
C[e]-lattice L CG. 
Output: The spectrum a = spec(L) G . 
(1) Compute Uq G GL^(C) such that 



\ 



\ 



Bqj 



where Bq G C^*^'^^ with spec(So) = i'^i} f^''^ ^ ^ [^^ 

(2) Set = ■= and M := U^^M. 

(3) Compute a minimal Grobner basis 

M := GB(M) G C[^]^^'' 
compatible with the ordering > on {O^cf)^ | G Z x defined by 

{k,i) > {l,j) ■.^k>ly {k^l M> j) 

for all {k,i), {l,j) G Z x 

(4 ) Return cr G Q*' with 

a{k + ai) := #({lead(M)}(fe,i)) 

for all {k,i) G Z x 
Lemma 18 Algorithm 3 terminates and is correct. 



Proof. Since —Td-,-(t> — (t>B, 

-Tdr{e^f) = 9''^{Bl + q) mod (^V)<C 

(k,i)<(i,j) 

with spec(SQ + I) — {aj + I}. Then, by Lemma 14, 

(i,fc)<(i,0 
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and hence, since M is a minimal Grobner basis, 

spec (L)(q;j + /) = dime S^^a+iiL/OL) 

= dime ((gr^ L/ 9 gr^ L)^^+^ 

= dime ((lead(M) ) C [6] /^(lead(M)) C [9] ) 

= dimc(((lead(M))C)(/,,)) 

= #({lead(M)}(,,,)) 

for all e Z X 

The following lemma reduces the problem of equality of C[^] -lattices to the 
problem of equality of filtrations on a finite vector space. 

Definition A G[6]-lattice L d G defines an increasing filtration L, on 
G by €>[9]-lattices Lp :— t^L and a corresponding decreasing filtration L* :— 
Ln^,. We denote the filtrations defined by Gq by G» and G* . 

Lemma 20 Let L G G be a C[9]-lattice. Then 
is an isomorphism for all a & Q and p & 'Z. 

Proof. This follows from 9PVaG = Va+pG and L"-^' = tPL for all a e Q and 
peZ. 

Definition 21 The sum J2 '■ IN*^ is defined by 

:= ^ aa{a) 

for allae IN^. 

The following lemma gives a criterion on the mean value of the spectrum to 
check equality of filtrations on a finite vector space. 

Lemma 22 Let F* and F* be decreasing filtrations on a finite vector space 
V with F* C F*. Then Z^spec(F2*) < Z^spec(F*) and equality implies that 

Proof. This is an elementary fact from linear algebra. 



15 



The following criterion on the mean value of the spectrum shall be used to 
check equahty of C[^]-lattices. 

Lemma 23 LetL2 G Li G G be C[9]-lattices. T/ien Espec(Li) < Espec(L2) 
and equality implies that Li — L2. 

Proof. Since L2 C Li, 

L;gr;o,i)GcL-gr|^,,)G 
where grj^ ^) G — 0o<a<i gr^ G and hence, by Lemma 22, 

Yl spec (l* gr[^ 1) G) <Y1 spec (l* gr J; 1) G) . 
By Lemma 20, spcc(Lj)(Q; + p) = spcc(^L* gr^ G^ {n — p) and hence 
spec{Li) = Y I] + ^ - spec (l* grj^ G) (p) 

0<a<lpGZ 

= n/x + ^ a dime (gr^ - ^ spec (l* gr^ G) 

0<a<l 

= n/x + ^ a dime (gr^ G) - ^ spec (l* gr g) 

0<a<l 

for i = 1, 2. This implies that 
^spec(L2) -^spec(Li) = ^ spec (l* gr ^ G) - ^ spec (l* grf^ 1) G). 

Let X e (Li\L2)n (^T4+pG\^y<a+j,G^ with < a < 1 and minimal a+p. Then, 

in particular, x ^ ^Li and hence, by Lemma 20, 7^ [t^x] e gr^^^gr^G. 

Moreover, there is a g > 1 such that G L2\9L2 and, again by Lemma 20, 
^ [rfx] e gr^^^'-^gr^G. This implies that L^'^gr^G C L^^gr^G and 
hence 

L-gr;^,,)GCL-gr|^,,)G. 
Then the claim follows from Lemma 22. 

The following theorem gives the mean value of the spectrum of Gq. 
Theorem 24 (C. Sabbah [1, 11.1]) ^Espec(Go) = 

By Theorem 24, one can compute t on Gq using Algorithm 1, 2, and 3 by 
increasing k until ^ X!spec(Lfc) = 

Our final goal is to compute a good basis of Gq. In terms of a good basis of 
Go, the matrix of t has degree one and its degree one part determines the 
spectrum of /. 
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Definition 25 Let (f) be a C[9]-basis of a t-invariant €[9]-lattice L C G. 
Then is called good if A± = A^+OAj where A^, Aj e C'^^'^, 



V 



a 



and (pi e VaiL for all i e [1, /i]. 

Lemma 26 Let (f) be a good basis of a t-invariant G[6]-lattice L <Z G and 



I 



0L\ 



\ 



\ 



a 



■.^A\. 



Then spec(L)(Q;) = e [1, //] | ccj = ol). 



Proof. Since is a C[^^]-basis of L and 0j G Vq,- for all i G [1, 

is a C-basis of gx^{L/9L) and hence spec(L)(Q;) = ^{i G [1, A*] | O-i = ol}. 



5 Monodromy 



Let Too be the monodromy of M around the discriminant D(/) = f{C{f )) of 
/ and To be the monodromy of G at 0. 

Theorem 27 (C. Sabbah [16, 1.10]) Too = f^^- 

Using Theorem 27, the monodromy T^o can be read off from the matrix of t 
with respect to a good basis. 

Proposition 28 Let 4> be a good basis of a €'[9]-lattice L G G and 



( 



\ 



a 



Then 



exp(-27ri(grr(4)+4 
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is a matrix ofToo where 

(K(C)),, := 

/or C = (qA.- e C'^x^ 



Ci,j if ai > aj + a, 
else, 



Proof. Since is a C[^^]-basis of L and (pi e V^^ for all i e [1,/^], 
is a C-basis of gr^{L/9L) and hence, by Lemma 20, 

-rdr = dtt : gii grr G ^ gr^ gr^ G © gr^"^ gr^ G 
for all a e Q and p e Z where ^A?" is induced by 

grf gr^, i = gr\ {^o o 0"^) = o gr];" (^5^) o 
Then, by [19, 6.0.1], exp(27ri(grf (^f) +Af)) 

is a matrix of Tq and hence, by 
Theorem 27, exp^— 27ri(gr]^^74^j + ^^)) is a matrix of T^q. 

6 Good lattices 

The following property is sufficient for the existence of a good basis of a C[^]- 

lattice [1, 5.2]. Recall that a morphism N : F'Vi >F'V2 of filtered vector 

spaces F'Vi for i = 1, 2 is called strict if Niyi) n F^ = N{F^) for all p e Z. 

Definition 29 We call a t-invariant C[9]-lattice L C G good if 

{rdr + a)P:L' gr^ G > L-p gr^ G 

is strict for all a & Q and p> 1. 

The following theorem follows from the fact that 

7V:= (rdr + a) 

0<a<l 

is a morphism of a natural mixed Hodge structure on the moderate nearby 
cycles = 0o<a<i g^^a ^ '^^^ Hodge filtration induced by G, as defined 

in Definition 19 [1, 13.1]. 

Theorem 30 (C. Sabbah [1, 13.3]) Go is a good C[9]-laUice. 
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The following lemma shall be used to construct an opposite filtration of L* on 
gr^ G for a good lattice L. 

Lemma 31 Let V be a finite vector space, F* a decreasing filtration on V 
with PP ^0 for p> m, and N e End(y) such that 

NP : F'V >F—PV 

strict for all p>l. Then Eg>o^^(^'") = 0<?>o^^(^'") and 

NP : F-(v/Eq>oN''{Fn) ^ F-p (v / Eq>o N'^ {F"^)) 

is strict for all p > 1. 

Proof. If X e F"^ with Np+\x) G Ifq=oN%F"') C F'^'P then 

NP+\x) e NP+\F"'+^) = 
since Np~^^ is strict and F"^~^^ — 0. Hence, 

7VP+i(F"^) + ^ Ar«(F'") = NP+^F"") (B^N^F"^) 

q=0 g=0 

and, by induction, Eg>o^''(^'") = ©9>o ^''(^'")- 

Let NP{x) e F« + Er>o N''{zr) with Zr e F"" for all r > 0. If m - p < g then 
NP(x - Y.r>pN''~^{zr)) e F'^-P+'^ and hence 

NP{x) e ArP(F'"+^) + 5] Ar'-(F'") C NP{F'i+P) ^Y.N'{F'^) 

r>p r>0 

since Np is strict and F"'+^ = 0. If m-p > q then A^^(x-Er>p ^''"^(-^r)) e F^ 
and hence iVP(a;) G Np{Fi+p) + Er>o N''{F"') since A^f is strict. This implies 
that NP is strict modulo Er>o ^''(-^'") for all p > 1. 

The following algorithm computes a C[t, ^]-basis of G compatible with the 
V-filtration refined by an opposite filtration of L' on gr^ G for a good lattice 
L. This basis shall be used to compute a good basis of L. 

Algorithm 4 

Input: (a) A matrix B = Ej>o-^i''"* ^ C[r]^^^ such that —rdrCt) = <j)B for a 
C[T]-basis 4> of and spec(i?o) = {a} with a > ai > ■ ■ ■ > > 
a — 1. 
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(b) A matrix M e C[t,9]^'^'' such that is a basis of a good C[9]- 
lattice L G G. 

Output: A matrix U = {U^''^){i,p)e[i,u]xz G GL^(C) such that {(f)U^''^)g>p is a C- 
basis of LPgTl.G and {{Ode - ai){^U'^P)} C + V^_i for all 

(1) Compute Uq G GL^(C) such that 



-^0 



V 



\ 



where G C^'^'^' with spec(i?Q) = for i e [1, u]. 
(2) Set = (0^)ig[i,^] := 0?7o and M := t/Q-^M. 
Compute a Grohner basis 

M := GB(M) e C[^]''^'' 

compatible with the ordering > on {9^(f)' | (p, i) G Z x [1,^^]} defined by 

{P, i) > {q,j) -.^ P > qy (p ^ q M > j) 

for all {p,t),{q,j) G Z x [l,z/]. 
(4) Set (Mf'^)(p,i)ezx[i,.] := M where 

{lexp(M^'0} = {(p,z)} 

for all {p,i) G Z X [1, 1/]. 
For i = 1. . . . .r do: 

(a) Compute {F'-'P)p^z e C*^''* 5?xc/i t/iai 

F^'P-.^ (T^lead(M«'0),<„-p. 

5etiV, -.= 31 -a,, 
(c) Compute U' = {U''P)pez e C^^^^^' stic/i i/ia^ 



/or all p G Z. 
Return 



U = (t/*'^)(i,p)g[i,,^]xz := ^7° 



V 



Lemma 32 Algorithm 4 terminates and is correct. 
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Proof. Since —rdrcf) — (f)B, 

-rdr{e''f) = e'^f{Bi + q) modV„_i 
with spec(SQ + q) — {aj + q} and hence, by Lemma 14, 

for all {q,j) G Z x [1, z^]. Since M is a Grobner basis, this implies that 0M^'-^ e 
Va^.+gL for all {q,j) G Z x [1, z/]. Then, by Lemma 20, 

L- grj:, G = {{fF^'-)C + K.)/V^<a, C gr^^ G 

and, since L is good and {Ode — C(i)4>^ = 4>'Ni mod Va-i, 

Nf : {F''')C >{F''—P)C 

is strict for all i e [1,^^] and p > 1. Hence, by Lemma 31, one can compute 
= ([/''f)pez e C/^^^"' such that 

(F^'f)C = ([/^''' \q>p)C, {NiW^P} C {[/^'f-^} 

for all e [1,1/] x Z. Then 

(^a^ - a,) i^W'^ = I'NiU^'P mod K_i 

and hence {{680 - ai){^U''P)} C {^U^'P-'^} + K-i for all {i,p) e [1, i/] x Z. 

7 Good bases 

The following algorithm computes a good basis of a good lattice L by a simul- 
taneous normal form computation and basis transformation. The computation 
requires a C[r, ^^]-basis of G compatible with the V-filtration refined by an op- 
posite filtration of L* on gr^ G. 

Algorithm 5 

Input: (a) A matrix B e C[r]^^'' with spec(i?o) = {«} and a > ai > ■ ■ ■ > 
> a — 1 such that —rdrcj) = (f>B for a C[T]-basis (f)ofVa 

(b) A matrix M e C[t,9]^'^'' such that 0M is a basis of a good C[9]- 
lattice L C G. 

(c) An indexing (p = (0*'^)(i,p)e[i,i/]xz such that (0*''^)g>p is a <C-basis of 
LPgrlG and{i9d^-a,)if'P)} C {f '^-'} + Kca, for all {t,p) G 
[l,u] X Z. 
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Output: A matrix M G C[t, 9]^^'^ such that (f)M is a good basis of L. 

(1) Compute a minimal Grobner basis 

M := GB(M) e GieY""^ 

compatible with the ordering > on {0^(j)''^ \ {k,i,p) G Z x x Z} 
defined by 

{k,i,p) > {l,j,q) -.^ k > I y {k = I A {i > j V {i ^ j Ap > q))) 

for all {k,i,p), {l,j,q) G Z x x Z. 

(2) Set (M*^'0(fe,,)6^x[i,.] M where 

{lexp(M'='^)} = {{k,i,n-k)} 



for all {k,i) G Z x [1, z/]. 
(3) Compute {A';'lj)(s,j,i)e^x[i,u]xz and 



- (a, + k) + ede)M''^ - M^+^'Mo^;i^,,, 

{V ,j')<{\+k,i) 
{s+l,j,n-l)<(l',j',n-l') 

such that lexp($^'[^^.) < (s + l,j,n - I) for all {k,i) G Z x for 
decreasing (s + l,j, n — I) until j = or ^^'/j 7^ and s > 1. 
(4) If^'lij = t/ien return M := (M'^'0(M)ezx[i,.]. 

("j; Choose {k, i) G Z x [1, i/] and {s,j, /) G Z x [1, z/] x Z wii/i 7^ and 
s > 1 suc/i that {s + n — I) is maximal. 

( 6) Set c'^'lj :— {1 + k + ai — s — I — Oij)"^ and 

(7) Go to (3). 

Lemma 33 Algorithm 5 terminates and is correct. 



Proof. By Lemma 14, 



for all {q,j) G Z x [1, i/]. Let m G {M)€[9] with lexp(m) = (A;, i,p). Then 

rVm e (f ''^ \q<p)C + V^^^G 
and, since 0m G {^M)€[9] = L and L"-'^ gr^. G = | g > n - A;)C, 

e (0^'^ I g > n - A;)C + y<aiG. 
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In particular, p > n — k. Moreover, T'^01ead(m) e gr^gr^. G and hence, by 
Lemma 20, 

^"-f-Vlead(m) e gr^gr^L = (01ead(M))C[^]. 

In particular, if p > n — k then lead(m) e 6'(lead(M))C[6']. Since M is a 
minimal Grobner basis, this implies that 

{lexp(M*^'^)} ^{{k,i,n- k)}, M'''' = lead(M'='^) mod terms < {k, i) 
for all (fc, i) e Zx [1, i/]. In particular, {0M*^''} C V^+ai for all (A;, i) e Zx [1, i/]. 
By Lemma 7, t o = 9[—Td-j) o = o 6'(5 + 6'c?0). Since 

{^(^9, - (A; + C {9'^^^^--^^+^)} mod y<i+fe+.„ 

there is a matrix ^o_i_|_fe,i ^^^^ ^^^^ 

9{B -{k + ai) + 9de)M^'' = M^+^^' A^^^^^ - mod terms < (1 + A;, i) 
and hence there are matrices A^sXj ^^^^ 

9{B -{k + a,) + 9de)M'^' - M^+^^'M^f,^,^, = ^ ^'M^''^'^'}, 

(s'+r,j')<(i+'=,») 

for all {k,i) e Z x [l,z/]. Choose {k,i) e Z x [1, z/] and (s, j, /) e Z x [1, z/] x Z 
such that (s + Z, J, n — Z) is maximal with ^s},j ^ s > 1. In particular, 
(s + /, j) < (1 + /c, i) and hence l-\-k-\-ai — s — I — aj >0 and c^'J^ > is 
defined. Moreover, since 

{{9de - (a, + C {^y>-(^+')} mod y<«.+,, 

= 0{B -{k + a,) + ^5,)(m^'^ - cy/^-iM''^A^;^^) 

{i',j')<{i+k,i) 

{s+l,j,n-l)<{l',j',n-l') 

= 9'M''^Al'l. + c';'l^9{9de -k- a,)0'-^M''^ A'^'I- 
= 9'M''^A^;l. + c]ij9'{9de + s-l-k- Q^i)M'^^ j^^^ 
^ 9^M'-^A^l + c£r (s + l + a,-l-k- a^M^^^A^i, 
= mod terms < {s + I, j, n — I) 

and hence {s + l,j,n — I) is strictly decreasing until ^'^'Jj — 0. Then the 
algorithm terminates and (f)M — (0M'^'*)(jfc^j)g^x[i,z/] is a C[^]-basis of L with 

ti^M''^) = IM'+''^A^'^^,^, + Yl + 0{k + a,)^M''^ 

{V ,3')<{l+k,i) 
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and {^M'^'^y C Vk+ai for all {k,i) e Z x [1, i/]. Hence, is a good basis of 
L. ~ 



The following algorithm combines Algorithms 1, 2, 3, 4, and 5 to compute a 
good basis of Go- 

Algorithm 6 

Input: A cohomologically tame polynomial f G C[x\. 
Output: (a) A vector (p e €[x,9]^ such that [0] is a good basis of Gq. 
(b) The matrix A = A^^ e €[9]"''^ oft with respect to [0]. 

(1) Setk-.^degif). 

(2) Compute G C[x, 9Y and A G €[9Y^^ by Algorithm 1. 

(3) Compute U G C[^]^''^ and B G C[r]^^^ by Algorithm 2. 

(4) Set := lU. B := U-\B - Tdr)U G €[tY, and M := [/"^ G C[r, 9^. 

(5) Compute a by Algorithm 3. 

(6) // ^ E 0- > ^ then set k :^ k + 1 and go to (2). 

(7) Compute U = (C/''^)(j,p)6[i,i/]xz £ GL^(C) by Algorithm 4- 

(8) Set I ~ (0C/''^)(^,p)6[i,.]xz, B := U-\B - Tdr)U G C[T]^x^ anc/ M := 
U-^M. 

(9) Compute M G €[T,9Y''^' by Algorithm 5. 

(10) Set := an(^ A := M-^9{B - Tdr)M G C[^]'^^'*. 
(^iij Return (p and A. 

Proposition 34 Algorithm 6 terminates and is correct. 



Proof. Let C Go be computed by Algorithm 1. Then = Go foi' ^ 3> 
and k is strictly increasing while > By Lemma 23 and Theorem 

24, L = Go if and only if ^E^x = ^Espec(L) = iEspec(Go) = This 
implies that = Go after finitely many steps. By Theorem 30, L :— — Go 
is a good lattice as required by algorithms 4 and 5. Hence, the algorithm 
terminates and is correct. 



Remark 35 In the local situation, one can replace the algorithms [8, 7.4-5] 
by the algorithms 4 and 5 to avoid the linear algebra computation [8, 7.4]. This 
modified algorithm is implemented in the Singular [20] library gmssing. lib 
[21]. 



24 



8 Examples 



Algorithm 6 is implemented in the Singular [20] library gmspoly.lib [22]. 
Using this implementation, we compute a good basis of Go for several exam- 
ples. By Lemma 26, the diagonal of Aj determines the spectrum of /. Using 
Proposition 28, we read off the monodromy Too around the discriminant of 
/ from A^. First, we compute two convenient and Newton non-degenerate 
examples [10]. 



Example 36 Let f ^ + + x^y^. Then Singular computes 



= (i,xy,y,x, 



and 



/_i i ^ ^ 

2 4 



-1 000 01000 
= -1 + 1 
000 -1 00010 
10 

The monodromy T^ has a 2 x 2 Jordan block with eigenvalue —1. 
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1 

"2/ 
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Example 37 Let f^x + y + z + x^y^z^. Then Singular computes 



e^x 



25 

Sex^ + a:^ ^x, m'^x^ - ^Ox-" + -^9x + —x' 



25, 

T 



25, 
T 



and 






25 
8 


"1 




^1 0^ 








125 
8 




10 


1 








+ e 


1 


1 










2 


1 









^0 1^ 



The monodromy T^ has a 2 x 2 Jordan block with eigenvalue 1 and a 3 x 3 
Jordan block with eigenvalue —1. 

Finally, we compute a non-convenient and Newton degenerate but tame [23, 
3] example. 
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Example 38 Let f = x{x^ + y^Y + x. Then Singular computes 



0= (l,623645y,^x,-2470(y=' + x2), 

- 11339(/ + x2y),475(2^y-^ - S^a;^ + Gx^), 3^V + V, 
6670 (^y^ - imx'^y + 6x^2/), S^^y^ - 20^^x^ - IS^x^ + 18x^ + 3, 

- 4365515(35^2?/^ - ?,m'^x'^y - ?,mdx^y + 180x^7/ + 24?/), — - — xy, 
- 8645(^x + 2y^ - 4x2), -124729(5^xy + - ^x^y) 



andA^t^^Af + eAf where 
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The monodromy T^ is unipotent with eigenvalues 



, „ , „ , „ ,1, 
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